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I. INTRODUCTION 

Opposite-signed two-muon production in deep inelastic 
scattering (DIS) with neutrinos or anti-neutrinos serves 
as a reliable signal for charm quark production through 
the sequence 

i^^N^^i-DX, D^fi+j^^X', (LI) 

where D represents a stable charmed hadron. This is 



especially true in dense targets such as the NuTeV detec- 
tor at Fermilab where absorption of pions and kaons in 
the hadronic shower suppresses backgrounds from meson 
decay that occur via 

i^^N ^ 11^71+ {K+) X; ^+ {K+) ^ (L2) 

Previous studies of DIS two-muon production in neu- 
trino interactions have yielded important measurements 
of the CKM matrix elements Vcd and Vcs, the effective 
charm quark mass, and the size and shape of the nucleon 
strange sea. These studies were performed in the context 
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FIG. 1. Schematic drawing of the Lab E detector. Beam 
enters from the left. The target/calorimeter is on the left and 
the toroid spectrometer is on the right. 

of the QCD-corrected quark-parton model and hence 
were restricted to large values of energy transfer, to 
the struck nucleon resulting in large amounts of observed 
hadronic energy, Ehad- 

At low hadronic energy, one expects to observe charm 
production via another mechanism, namely, diffractive 
production of pseudoscalar and vector Ds mesons via 
both coherent and incoherent scattering: 



1^ 



Vi^A 



(1.3) 



where the low energy decay photon 7soft will accom- 
pany production of vector D*g, A = Fe for coherent pro- 
duction in our experiment, and A — n or p for incoherent 
production. Here, large energy transfer to the meson may 
not result in large hadronic energy because the nucleus 
remains intact. 

These reactions have been observed at the few event 
level in bubble chamber and emulsion experiments . 
It is important to understand the size of this diffractive 
contribution because of its influence as a background to 
DIS charm production. These processes are also of in- 
terest to future high statistics neutrino experiments at a 
muon coUider. They provide the possibility of measur- 
ing the ratio Vcd/Vcs via comparison of the rates Vf^A 
fi^D*+A and v^A ^"Dp'A. AdditionaUy, they can 
create a signature which mimics quasi-elastic production 
of T— leptons through the chain i^f^A /i"!)^ (tsoft) A, 
t^Vt , which may be of concern to high sensitivity 



Vr oscillation searches. 



A competing reaction that produces the same experi- 
mental signature, an opposite-signed muon pair with van- 
ishing Ehad , is the neutrino trident process P,B| : 



,A- 



A. 



(1.4) 



In principle, the neutrino trident reaction provides an in- 
teresting test of electroweak theory since contributions 
from W and Z decay produce a reliably calculable 40% 
destructive interference effect. In practice, the very small 
cross-section implies that only a handful of neutrino tri- 
dents have previously been observed in neutrino scatter- 
ing. Furthermore, the neutrino trident process must be 
considered in combination with the expected signal from 
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FIG. 2. Drawing of one segment of the Lab E calorimeter. 

diffractive charm production in experiments that are only 
sensitive to two-muon final states. This point has not 
been recognized in previous measurements of neutrino 
tridents. 

NuTeV (Fermilab Experiment E815) is a high- 
statistics experiment studying neutrino and anti-neutrino 
interactions with a high intensity, high energy sign- 
selected beam. Its primary goal is to measure the 
electroweak mixing angle sin^^vy; however, the high- 
statistics nature makes it possible to study the rarer pro- 
cesses described above. At low (zero) hadronic energy the 
dominant deep- inelastic scattering (DIS) processes fall off 
and additional channels such as diffractive / D*^ pro- 
duction H] , neutrino trident , and diffractive J/ip ^ 
become more important. This is the first analysis to 
study the inclusive production of all processes in the final 
state with low hadronic energy and two muons. 

The next section will describe the NuTeV experiment 
while Section III details the data selection. Predictions 



for the sources of Xovj-Ehad two muons events are given 
in Section [V. Section ^ details the charged- and neutral- 
current analyses and conclusions are stated in the final 
section. 



II. THE NUTEV EXPERIMENT 



The NuTeV experiment ran at Fermilab using the re- 
furbished Lab E detector |^ and and i/^j beams pro- 
vided by the newly installed Sign-Selected Quadrupole 
Train (SSQT) §. The SSQT has the capability of select- 
ing either muon neutrino or muon anti-neutrino beams 
while leaving ^ 2 x 10~^ of the anti-selected type. 

The Lab E detector, located 1.5 km downstream of the 
primary target, consists of a target/calorimeter followed 
by a toroid spectrometer (Fig. M. The calorimeter is 
composed of 42 segments of four steel plates, two liquid 
scintillator counters (SC) and one drift chamber (DC) 
(Fig. H). The steel provides mass for the neutrino target, 
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the SC measure the longitudinal vertex and energy de- 
position {E), and the DC are used to find the transverse 
vertex 

and reconstruct downstream tracks. Because of the 
high density of the target, only muons deeply penetrate 
the calorimeter, and all other particles create a hadronic 
or electromagnetic shower near the interaction vertex. 

The toroid spectrometer, located immediately down- 
stream of the calorimeter, focuses muons from the pri- 
mary charged-current vertex given the type of beam (z^^ 
or D^), measuring both charge and momentum (p) of 
muons with p > 5 GeV/c which enter the spectrome- 
ter. It is also possible to measure the momentum of a 
subset of muons with 5 < p < 15 GeV/c which range out 
in the calorimeter, and to place a lower bound on the mo- 
mentum of muons which exit the side of the calorimeter. 

The NuTeV detector was calibrated with a separate 
beam of hadrons, muons, or electrons throughout the 
running period of the experiment. Hadronic and muon 
energy scales for the calorimeter were determined to 
0.43% and 1.0% respectively over the energy range 5—200 
GeV 1^. Muon momentum measurement in the spec- 
trometer is limited by multiple Coulomb scattering to 
Ap/p = 0.11, and the sampling-dominated hadronic res- 
olution was approximately AE/E = 0.86/^/E. 



• The second muon's momentum could be obtained 
by either toroid spectrometer or range information; 
its momentum was required to be at least 5 GeV/c. 

• Either the calorimeter drift chamber or scintillator 
counter signals had to be consistent with passage 
of two muons over five drift chambers or seven scin- 
tillator counters starting at the event vertex. 

IV. LOW HADRONIC ENERGY, TWO MUON 
SOURCES 

The analysis strategy consists of comparing data to a 
model comprised of all known sources of events with two 
muons and small hadronic energy. Kinematic distribu- 
tions were generated according to electroweak theory for 
neutrino trident production, leading order quark-parton 
model predictions for DIS feed-down, and vector meson 
dominance (VMD) or partially conserved axial current 
(PCAC) models extended to four flavors for diffractive 
charm production. Detector response was modeled using 
a GEANT-based Monte Carlo (GEANT 3.21), with sim- 
ulated events processed in an identical fashion as data. 
The following sections describe the various sources and 
predictions for their rates. 



III. EVENT SELECTION 



A. DIS Charm Production 



Analyses presented here use the full NuTeV data sam- 
ple from the 1996-1997 flxed-target run corresponding to 
1.3 X 10^^ protons on target (POT) for neutrino running 
and 1.6 x 10^^ protons on target for anti-neutrino run- 
ning. The size of various event samples from NuTeV are 
listed in Table |. 

Events were selected for the Iow-Ehj^^jj two-muon anal- 
ysis based on the following criteria: 

• The event vertex was required to be at least 25 
cm from any side, 40 cm steel-equivalent from the 
upstream end and 1.4 m steel-equivalent from the 
downstream end of the calorimeter. 

• At least two muons were required to be found and 
fitted by the tracking code. At least one of these 
had to be reconstructed in the toroid spectrometer 
with a momentum greater than 9 GeV/c. 



TABLE I. NuTeV event samples. 





# of events 


^ of events 




(^) 


(^) 


Single muon (all Eh ad) 


1.3 X 10" 


0.46 X 10" 


Two muons (all Eh ad) 


4300 


1300 


Single muon {Eh ad < 3 GeV) 


0.10 X 10^ 


0.06 X lO'' 


Two muons {Ehad < 3 GeV) 


33 


21 



Neutrino DIS produces two major sources of two-muon 
events. Charged-current charm production can result in 
a second muon from the decay of the charmed meson 
(Fig. 11(a)). Also, during charged-current interactions, a 
t:/K in the hadronic shower can decay to a muon and a 
neutrino before interacting (Fig. §(b)). 

Two-muon events from DIS charm were modeled via 
leading order predictions using d and s parton distribu- 
tion functions measured in this and previous experiments 
P-p^ assuming an effective charm quark mass rric = 1.32 
GeV/c^. Charm quark fragmentation was treated using 




(a) 



(b) 



FIG. 3. The Feynman diagrams for DIS production of 
two-muon events: (a) DIS charged-current production of 
charm with a semi-muonic decay; (b) DIS charged-current 
production with a n/K decay in the hadronic shower. 
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FIG. 4. Distribution of expected hadronic energy (Ehad) 
for DIS production of two-muon events (Monte Carlo). This 
analysis is concerned with the contribution near zero. 

the CoUins-Spiller scheme jl^ with fragmentation pa- 
rameter e — 0.93 and : : A^t production ratios 
measured in Fcrmilab E531 and corrected for effects 
of Dg production |l6| . Charmed hadron decay was mod- 
eled using the best available data as summarized by the 
Particle Data Group , with particular attention given 
to purely leptonic decays. 

Two muon events from n/K decay in the hadron show- 
ers of ordinary charged-current (CC) events were mod- 
eled using measured parameterizations from a previ- 
ous experiment (FNAL E744/770) which ran at similar 
neutrino energies and used the same neutrino target. 

The expected hadronic energy (Ehad) spectrum for 
DIS two-muon Monte Carlo events is shown in Fig. ^. 
While these sources generally contain significant amounts 
of hadronic energy, some contribution at low E^ad is 
seen. Normalizing to the high-En ad (Ehad > 5 GeV) 
two-muon sample we predict 10.4(4.0) events with Ehad 
< 3 GeV in ^(i') mode due to the low hadronic energy 
tail of the DIS process. 




(b) 

FIG. 5. The Feynman diagrams for neutrino trident pro- 
duction: (a) charged-current production; (b) neutral-current 
production. 




(a) (b) 

FIG. 6. The Feynman diagrams for (a) and (b) Dg 
production. 



B. Neutrino Tridents 

Neutrino trident production is a purely electroweak 
process in which interference between the charged- and 
neutral-current diagrams causes a 40% decrease (from V~ 
A) in the total cross-section |^,^,|l^ . Feynman diagrams 
are shown in Fig. ||. The small (~ 10"'* fb/nucleon) 
cross-section has limited the observation until recent 
years 

For this analysis, the Monte Carlo was generated us- 
ing the full matrix element with W-Z interference [ p3|j24| 
including contributions from both coherent nuclear and 
incoherent nucleon scattering. This procedure incor- 
porated all possible kinematic correlations between the 
two muons and represents an improvement over previous 
methods |20| , ^ . In particular, there is a strong correla- 
tion between the energies of the two muons which is very 
important to the acceptance: when one muon's momen- 
tum is high, the other is preferentially very low. 



The lengthy expression for the cross-section can be 
found in the references |,|,|l|,|3|j24|. The rate de- 
pends on the electroweak mixing angle, which we set to 
sin^ 9w = 0.2222, and, weakly, on nuclear and nucleon 
form factors. Standard dipole parameterizations were 
used for the latter with a vector pole mass rriy = 0.71 
GeV^/c'*. The nuclear form factor for iron was calcu- 
lated assuming a Fermi charge density function for iron 
with nuclear size parameter c = 3.9 fm and thickness 
parameter b = 0.55 fm. 

NuTeV should observe 4.8(2.2) neutrino trident events 
in J^(S') mode from both coherent and incoherent produc- 
tion. 



C. -D|/Ds Production 



Diffractive production 




has 



been observed in previous experiments [pjO] and con- 
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tributes to the low-E^f^^i two-muon sample whenever 
the Ds meson decays either to all-lepton final states 

or to final 



states with small hadronic energy (_Dj 
The latter case was modeled by assuming that the L'J 
semi-muonic decay rate was saturated by the channels 



in the measured propor 



tions summarized by the Particle Data Group ||17|. Only 
contributions from coherent diffractive production were 
considered. Incoherent production of D*^) from vN 
scattering is already included in the inclusive DIS rate. 

Vector D*g production was modeled assuming a VMD 
type mechanism with cross-section given by 



d^aiy^N -> fiD*N) 



Ml, 



{l-e) 



(4.1) 



where and v are the momentum and energy transfer 
to the D*g , t is the square of the momentum transfer 
to the nucleus, gp is the p coupling constant, M^j is 
the mass of the D*g meson, E is the incoming neutrino 
energy, e is virtual W polarization (e = {AE{E — v) — 
Q^)/{4:E{E - i^) + Q2 + 2i/2)), and b is the slope of the 
distribution of momentum transfer squared (t) to the nu- 
cleus (nucleon) (6 = 3 incoherent; b — 145 coherent). A 
check was performed which showed the NuTeV experi- 
ment is insensitive to effects of the virtual- polarization 
on the Dg decay. The overall normalization was left to 
be determined by the data. To set the scale, the normal- 
ization was also obtained in the SU (4)-flavor limit by 
comparing to the measured cross-section for diffractive 
p"*" production in neutrino scattering. 

The number of events expected by NuTeV is normal- 
ized to the observed inclusive two-muon data with Eh ad 
> 5 GeV 

^ X aM{E,) X e{E,) x BFm 
/ X af,f,iE,) X eiE,) x BF^,f, d<P 

(4.2) 

where M refers to the meson being produced, /i/i refers 
to DIS two-muon data, $ is the NuTeV flux, a is the 
cross-section, Ei^ is the incident neutrino energy, e is the 
detection efficiency, BF is the (semi-)muonic branching 
fraction, and A^^^ is the number of DIS two-muon events 
observed in the NuTeV data. This procedure predicts 
33.0(13.6) observed Dg events for i/(P) mode in NuTeV. 

The contribution of pseudoscalar Dg decay was esti- 
mated using PCAC formulas |2qJ2^ adapted for charm: 



/if 



dQ^dvdt 



^2 (l-e) (Q' 



Ml Y 




FIG. 7. The Ehad distributions for (a) D% and (b) D*g^ 
Monte Carlo. 

where M^g^ is the mass of the E)s\ meson and /ug is the 
Ds decay constant (/ds = 0.31 GeV). 102.9(48.2) events 
are predicted for v(y) mode in the SU(4) flavor limit. 

Figure |^ shows the Eh ad distributions predicted by 
the Monte Carlo. The r decay mode results in the peak 
at Ejhad = for production while the extra photon 
from the D*^ decay washes the peak out. 



D. Other Sources 

Several other contributions to Xov^-Ehad two-muon 
states contribute at the sub-event level. These include 
decays from neutral-current (NC) diffractively pro- 
duced J/V', V^'^ decays from NC diffractively produced 
vector mesons, CC diffractively produced tt"*" in which the 
pion decays, tt"*" from CC baryon resonance production 
followed by pion decay, and quasi-elastic CC scattering 
where pattern recognition errors split the outgoing muon 
track into two. 



1. NC Diffractive Vector Meson Sources 

Diffractive production of vector mesons is an impor- 
tant process in photoproduction experiments. A similar 
process is available in the weak sector with the substi- 
tution of a for the photon (Fig. ||). Vector mesons 
which can decay to two muons include , uj, (j> and J/ip- 




(4.3) 



FIG. 8. The Feynman diagram for neutral-current diffrac- 
tive production of light vector mesons (V°). Here V° can be 
a, , Lo, (j) or J/ip meson and N represents either a nucleus 
(coherent) or nucleon (incoherent). 
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The VMD model is used to predict the expected num- 
ber of vector mesons. Cross-sections are normahzed 
to the Fermilab E632 measurement of neutrino produc- 
tion of single |^ with the neutral-current (NC) and 
charged-current (CC) cross-sections related by Pumplin's 
procedure . The J/V' estimate is given in the SU (4) 
flavor limit for comparison purposes. 

Table || lists the number of vector meson events pre- 
dicted to be observed in the NuTeV data. All sources ex- 
cept J I'll) production are expected to be small, with the 
suppression largely attributable to the very smaUji^/i^ 
branching fractions for vector mesons. Section VB de- 



scribes a measurement of the J /ip signal from diffractive 
production. 




FIG. 9. The Feynman diagrams for single vr* production: 
(a) diffractive and (b) resonance. 



2. Sources from Pion Decay or Mis-reconstruction 

Single pious produced in NuTeV have two primary 
sources: diffractive and resonance production (Fig. ||). 
The density of the NuTeV detector limits the contribu- 
tion of both sources to the two-muon sample since the 
pion must decay prior to interacting. Diffractive produc- 
tion is predicted using the PCAC model normalized to 
the measured cross-section resulting in 0.08(0.02) es- 
timated events in i^(i^) mode. Resonance production can 
result in a decay of a final state pion with little visible 
hadronic energy. This is calculated with a model from 
Rein and Sehgal and predicts < 0.1 events in either 
mode. 

Quasi-elastic production of events (Fig. |l^) has 
also been considered ||3^. The muon from the A^ de- 
cay is of very low momentum resulting in extremely low 
acceptance. We predict <0.1 events in either mode. 

The final source considered was leakage from low- 
^HAD single muon events due to mis-reconstruction. 
Contributions from this class of events was greatly re- 
duced by requiring that the scintillator counters or 
drift chambers were consistent with two muons over 
a minimum length. The primary remaining source of 
these events was quasi-elastic scattering where multi- 
ple coulomb scattering and drift chamber inefficiencies 
caused the pattern recognition software to find a spuri- 
ous second track in the event. This contribution was es- 
timated using straight-through muons which result from 




FIG. 10. The Feynman diagram for diffractive A,-, produc- 
tion. 

charged-current interactions upstream of the detector 
that enter the calorimeter. A random longitudinal event 
vertex was chosen, and all event information upstream of 
this point was discarded, leaving a topology identical to a 
quasi-elastic event. These straight-through muons were 
analyzed in the same way as the neutrino data sample 
and the results normalized to the observed number of 
single muon quasi-elastic events. The resulting predic- 
tions are 0.75(0.25) events in i/(i/) mode. 

A summary of the predicted number of events from 
each source is shown in Table [II. The remainder of this 



paper describes the measurement of the largest sources: 
DIS, neutrino tridents, / D^^ , and J/4'- 



V. RESULTS 



TABLE II. VMD predictions for the expected number of 
diffractively produced light vector mesons which decay to two 
muons. 





# of events 


# of events 


^ of events 


^ of events 




(coherent) 


(incoherent) 


(coherent) 


(incoherent) 




(^) 


(^) 


(^) 


(^) 




0.0074 


0.0039 


0.0037 


0.0019 




0.023 


0.016 


0.011 


0.008 




0.036 


0.022 


0.018 


0.011 


J/4 


4.88 


17.1 


1.86 


7.76 



A. Charged-Current Analysis 

Our charged-current analysis examines three kinematic 
variables: the hadronic energy {Eh ad) up to 15 GeV, 
the two-muon invariant mass (M^^), and the absolute 
value of the smallest difference between the muon az- 
imuthal angles (A0). Distributions for these quantities 
as measured in NuTeV data are shown in Figs. |l|, |l^, 
O, and 14, Since the analysis is primarily concerned 
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TABLE III. Predictions for the number of events expected to be observed by 
NuTcV in the \ow-Ehad (Eh ad < 3 GeV) two-muon sample. 



# of events 


# of events 







DIS 

Neutrino Tridents 

^1 
Df 

Light vector mesons (p'', uj, (f>) 
J/ip (coherent) 
J/i/i (incoherent) 
Single TT* 

Mis-identified single muon events 



Sum 



10.4 
4.8 
102.9 
33.0 
< 0.12 
4.88 
17.1 
<0.1 
<0.1 
0.75 
174.4 



4.0 
2.2 
48.2 
13.6 
< 0.08 
1.86 
7.76 
<0.1 
<0.1 
0.25 
78.2 



RefTj:! 
Ref. |:4| 
Ref. |T "l 
Ref. Jl i 
Ref. M 2^ 
Ref. [26 ~ 
Ref. [26 25 
Ref. IIP 
Ref![BOr 
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FIG. 11. The hadronic energy distribution for the 
two-muon data. The first twelve bins contribute to figures |l3| 
and0. 




3 12 3 

A<p (radians) &<p (radians) 

FIG. 12. Distributions for the two-muon NuTeV events 
with Eh AD < 15 GeV. The plots are the two-muon invariant 
mass (Mfifj) for v (a) and v (b) modes and At/) in v (c) and v 
(d) modes. 
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FIG. 13. The two-muon invariant mass (M^i^) for the low- 
Ehad [Eh ad < 3 GeV) two-muon data. The curve shows a 
fit to the distribution x" e^'^+''"''> . 
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high Eh AD, and at A4> = and tt, the latter feature re- 
flecting the essentially two-body final state of DIS charm 
production. The neutrino trident signal peaks at low 
Eh AD, low Mfj_^ and is more uniform in A(f>. The Dg 
and Dg distributions are intermediate between DIS and 
neutrino tridents, with Ds distinguishable from Dg in 
the Eh AD distribution dvie to the decay photon contri- 
bution to Eh AD for D*g. Diffractive Dg and neutrino 
trident distributions are very similar, hence their mea- 
surements are highly correlated. 

Data and Monte Carlo were binned 
three-dimensionally in {Ehad, M^^, A4>) space: 18 bins 
in Ehad (-3 - 15 GeV), 6 in M^^ (0 - 6 GeV); and 6 
in |A(/)| (0 - TT radians). The Monte Carlo is fit to the 
data using a maximum-likelihood technique. The Monte 
Carlo sets are summed together using a single normal- 
ization factor for each source which was defined to be 1.0 
for the level presented in the preceding section. Neutrino 
and anti-neutrino modes were fit simultaneously. 

The results of the fit are shown in Table IV . Figure |l^ 
compares these results to the data and shows we are 
able to describe the data with the four largest sources. 
The DIS contribution is consistent with that expected 
from the higher Yjhad two-muon analysis The neu- 
trino trident contribution is consistent with the Standard 
Model prediction, but can not distinguish between V-A 
and the Standard Model. This is in contrast to previ- 
ous analyses which ruled out V~A but did not consider 
diffractive sources. 



FIG. 14. The Aiji distributions for the Iow-Ehad (Ehad 
< 3 GeV) two-muon data. At/i is defined in the text. 



with events that peak at low Ehad, the Af^^ (Fig. p^ ) 
and A(/) (Fig. [ij) distributions are plotted with the ad- 
ditional restriction Ehad ^ 3 GeV to emphasize the 
diffractive region. However, all events with Ehad < 15 
GeV are used in the subsequent fitting. Thirty-three 
neutrino and 21 anti- neutrino two-muon events were ob- 
served with ^HAD < 3 GeV. These numbers are signifi- 
cantly below the sum of the predicted number of events 
in Table III, indicating the naive SU{A) flavor predictions 
for heavy quark final states are too large. On the other 
hand, there is a large excess over the predictions from 
DIS feed-down and neutrino trident production. 

A fit was performed allowing up to four sources of low- 
E_ffAD two-muon events: DIS, neutrino tridents, diffrac- 



tive of, and diffractive Dg^. These sources were sim- 
ulated via Monte Carlo and normalized to the NuTeV 
inclusive two-muon samples. Figure |l^ shows distribu- 
tions of Ehad, M^^, and At/) for the four major sources. 
Deep-inelastic scattering contributions peak at high M^^, 



TABLE IV. Parameters from the three parameter fit to the 
Iow-Ehad two-muon data. 



Parameter 



Resuh 



DIS 

Neutrino Tridents 

Diffractive Charm (D| + £)J±) 



+0.09 
0-90 _o.08 
+1.73 
-0.72 
+0.06 
0-18 -0.06 



0.72 
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FIG. 15. Monte Carlo distributions of visible hadronic energy, two muon invariant mass and A(p for the four 
the four largest sources of low-E_frAi3 two- muon events: (a) DIS charm; (b) neutrino tridents; (c) diffractive D^; 
(d) diffractive D^^ . The invariant mass and A0 distributions are for Eh ad < 3 GeV. 



TABLE V. Measured event contributions to the Io\y-Ehad two-muon event sample. 

# of events # of events 

M M 



J/ip (90% CL) 


<7.5 


<5.0 




+0.9 


+0.4 


DIS 


9-4 -0.9 


3.6 _o.4 




+8.3 


+3.8 


space0.2cm Neutrino Tridents 


3-5 -3.5 


1-6 -1.6 




+8.1 


+3.7 


space0.2cm (D| + D^^) 


24.5 _8.2 


11 1 -3.7 



vspace0.2cm Mis-identified single muon events 0.75 ± 0.43 0.25 ± 0.25 
vspace0.2cm All other sources (estimate) <0.1 <0.1 
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FIG. 16. Comparison of the final result (MC) to the 
Iow-Ehac two-muon data for (a,b) F,had, (c,d) M^+^-, (e,f) 
A(/!>. The left side is u mode; the right side is u mode. 
The M^+^- and A<j) distributions are for Ehad < 3 GeV. 
The points represent the data while the histogram shows the 
Monte Carlo. 



The consideration of all sources of \ow-Eihad two- 
muon events allows us to measure diffractive charm pro- 
duction. The and D*^ sources have been combined 
in proportion to the theoretical predictions and a single 
fit parameter used. This yields cross-sections of 



(7 



^l^{Ds + D*s)Fe) =(3.3±1.1) fb/nucleon, 



evaluated at E^, — 130 GeV using the modified 
VMD and PCAC predictions to extrapolate in en- 
ergy under the assumptions a {uf^Fe 



-D*+Fe) 



a 



fi+D%-Fe) 



-D+Fe) 



and a {vf^Fe — > /i 
CT [pfj^Fe ~^ fi~^Dg Fe) . A second fit performed with 
the neutrino trident parameter fixed to the Stan- 
dard Model prediction yielded the consistent results 
a{Vf,Fe^ H'iDs + D*s)Fe) = (3.0 ± 0.8) fb/nucleon 
at E,y = 130 GeV. The quoted errors are completely dom- 
inated by statistics. This result assumes an isotropic 
Dg decay. Studies showed effects of a possible Dg po- 
larization to be small. The largest change, correspond- 
ing to nearly complete longitudinal polarization, lowered 
a{Ds + D*g) by 0.4 fb/nucleon. 

Previously, the Big Bubble Chamber Neutrino Collab- 
oration combined various data samples to measure the 
diffractive rate of charmed strange mesons ( + D*g^) 
per charged-current vl (/ is an isoscalar target) interac- 
tion g. They measured a rate of (2.8 ± 1.1) x 10"^. 
The observation of Dg^ production by CHORUS is in 
agreement with this rate. Using the results of our second 
fit, we find a rate of (3.2 ± 0.6) x lO'^, which is 



250 - 



200 - 



150 - 



100 - 



50 - 




12 3 4 5 6 
(GeV/c') 

FIG. 17. The two muon invariant mass (A/^^) for the J/ip 
Monte Carlo. The curve shows a Gaussian fit. 

consistent with previous results. 

Table ^ lists the number of events contribution of each 
source in the \ow-Ehad two muon data sample as deter- 
mined by this analysis. 



B. Neutral-Current Analysis 

Neutral-current J/ip production produces a clear sig- 
nature in the two muon invariant mass, particularly if 
Ehad < 3 GeV is imposed to select diffractively pro- 
duced events. There is no evidence for a J/tp signal in 
Fig. ^sj; however, the relatively poor resolution of the 
NuTeV detector may be obscuring a contribution from 
this source. To assess this possibility, a diffractive J/tp 
sample was simulated via Monte Carlo to obtain the M^^ 
distribution shown in Fig. ^ A Gaussian fit to this dis- 
tribution yields a resolution (Jq = 0.40 GeV/c^. 

A maximum likelihood fit was then performed to de- 
termine the amount of J/ip present in the data. The fit 
function was taken to be 



(5.1) 



where Af^^ is the two muon invariant mass. Mq and do 
are the mass and width of the J/^ as measured by the 
Monte Carlo. The first term represents a smooth param- 
eterization of the background description where a and 
7 determine the shape and (3 the normalization. The 
second term is a Gaussian description of the J/tp con- 
tribution with mean mass Mq and width tJo set to the 
Monte Carlo prediction. The parameter A measures the 
amount of J/ip in the data. 

The results of the fit are shown in Tabic 0. A 90% 
confidence level (CL) on the J/tp contribution is set by 
fixing the J/ip amplitude to various increasing levels 
and fitting for the background. The likelihood function 
(£(A)) was plotted as a function of A and the 90% CL 

limit 



set by J^^'"- £{A) dA/ C{A) dA = 0.90. The 
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FIG. 18. 90% confidence level limit on diffractive J /tj) pro- 
duction. The curve shows the fit to the background with the 
90% CL contribution from J /^ii. 

resulting number of events is found by integrating the 
Gaussian with amphtudc Acl and converting to a cross- 
section by normahzing to the DIS two-muon sample. The 
90% confidence level limits on the diffractive J/i/" cross- 
section are 

a {Vf,Fe v^J/tpFe) < 0.21 fb/nucleon at 90% CL, 

a {D^Fe D^J/ipFe) < 0.36 fb/nucleon at 90% CL, 

at a mean production energy of E^, — 175 GeV for 
Eh AD < 3 GeV. Figure |lj shows the results of this limit 
including the J/ip contribution. 

Using the VMD model, we can extrapolate the J/ip 
cross-section to lower E^, to compare to the CDHS mea- 
surement of (0.042 ± 0.015) fb/nucleon. Figure |l9| shows 
a 90% CL limit for the J/tp cross-section normalized to 
the NuTeV measurement and compares to the CDHS re- 
sult. The energy dependence of this limit is dependent 
upon the model inputs. We show the limit for a momen- 
tum transfer squared (t) distribution of e"^''^* j2^ and 
coherent production only. The interpretation of the 




50 100 150 200 250 300 

E„ (GeV) 

FIG. 19. Energy dependence of the 90% confidence level 
limit on diffractive J/tp production normalized to the NuTeV 
measurement. The solid line is the limit set with Ehad < 
3 GeV and the dashed line with Ehad < 10 GeV. The data 
point shows the CDHS measurement. 

CDHS measurement as diffractive production is contra- 
dicted by our limit; however, the CDHS Ehad cut of 
10 GeV compared to a mean neutrino energy of 70 GeV 
could have accepted DIS production of J /ip- 

The analysis was repeated with a higher ^had cut 
(^HAD < 10 GeV) which matches the CDHS selection. 
The results are also included in Table VL The limit as a 



function of E^, is shown as the dashed curve in Fig. |19| . 
While this limit is near the CDHS measurement, the 
diffractive J/ip Monte Carlo shows the lower Ehad (<3 
GeV) is appropriate for coherent diffractive production. 
One possible explanation is that the higher cut allows a 
contribution from DIS J/ijj production. 



TABLE VI. Limits on diffractive J/tp production observed by NuTeV in the 
Iow-Ehad two-muon sample. 







v mode 


1/ mode 


Ehad < 3 GeV 










# of events (fit) 


3.3 ± 5.0 


-1.7 ± 4.0 




# of events (90% CL) 


7.5 


5.0 




Average E^ {J/iP) (GeV) 


185.0 


168.0 


Cross-Section (90% CL) (£^ 


= 175 GeV) (fb/nucleon) 


0.19 


0.32 


Cross-Section (90% CL) {E, 


. = 70 GeV) (fb/nucleon) 


0.011 


0.017 


Ehad < 10 GeV 










# of events (fit) 


7.8 ± 14.3 


-2.4 ± 10.5 




# of events (90% CL) 


24.8 


11.5 


Cross-Section (90% CL) (£^ 


= 175 GeV) (fb/nucleon) 


0.63 


0.73 


Cross-Section (90% CL) {E, 


. = 70 GeV) (fb/nucleon) 


0.034 


0.040 



CDHS result (Ehad < 10 GeV) 

Cross-Section (Fit) {{E^ = 70 GeV) (fb/nucleon) (0.042 ± 0.015) 
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FIG. 20. The two muon invariant mass distribution for 
high-E/jAD (Eh AD > 10 GeV). The curve shows the fit to 
the background with the 90% CL contribution from J/ip. 

For completeness, the Ehad cut was changed to 
Eh AD > 10 GeV to allow a search for inclusive DIS J/ip 
production. Figure ^ shows the resulting Af^^ distribu- 
tion. A fit was performed with the background modelled 
by an asymmetric Gaussian function and the data binned 
in 0.5 GeV/c^ intervals. No statistically significant J/ij; 
signal was found, and 90% confidence level limits for the 
(DIS J/ip)/{DlS CC charm) rates were found to be 



< 0.024 at 90% CL, 



< 0.069 at 90% CL, 



when averaged over the NuTeV beam spectra. Assuming 
a similar energy dependence to DIS CC charm produc- 
tion, this corresponds to cross-section limits evaluated at 
= 125 GeV, 

cr {i^f,N I'^.J/ipX) < 2.2 fb/nucleon at 90% CL, 



cr {D^N ^ v^J/tpX) < 3.4 fb/nucleon at 90% CL, 

at an average DIS CC charm production energy of i?^ 
= 125 GeV. Figure |o| shows the 90% CL curve for the 
inclusive ( Ehad > 10 GeV) DIS two-muon sample. 



VI. CONCLUSIONS 

We have performed an inclusive analysis of the low- 
^HAD two-muon data sample available from NuTeV. All 
known Standard Model processes have been considered, 
and four significant sources were found to contribute: 
DIS, neutrino trident, diffractive Dg, and diffractive 



D 



*± 



We have measured diffractive Da and D 



*± 



cross-sections to be 1.4 ± 0.4 fb/nucleon and 1.6 ± 0.5 
fb/nucleon, respectively, for i/^Fe scattering at E„ = 130 
GeV, in agreement with previous results assuming a 



VMD energy dependence of the cross-section. We ob- 
serve a statistically weak neutrino trident signal consis- 
tent with Standard Model predictions. Finally, we see no 
evidence for either diffractive or DIS production of J/^p 
in NC and scattering on iron. 
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